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Abstract: We present a complete study of the Ergodic Theorem for the difficult prob- 
lem of Non-Linear Wave Propagations through Cylindrical Measures/Path Integrals and 
the famous Ruelle-Amrein-Geogerscu-Enss (R.A.G.E.) theorem on the caracterization of 
continuous spectrum of self-adjoint operators 

1 Introduction 

One of the most important phenomenon in numerical studies of the non-linear wave motion 
especially in the two - dimensional case - is the existence of overwhelming majority of wave 
motions that wandering over all possible system's phase space and, given enough time, 
coming as close as desired (but not entirely coinciding) to any given initial condition ([1]). 

This phenomenon is signaling certainly that the famous recurrence theorem of Poincare 
is true for infinite dimensional continuum mechanical systems like that one represented 
by a bounded domain when subject to non-linear vibrations. 

A fundamental question appears in the context of this recurrence phenomenon and 
concerned to the existence of the time average for the associated non-linear wave motion 
and naturally leading to the concept of a infinite - dimensional invariant measure for the 
non-linear wave equation ([1]), the mathematical phenomenon subjacente to the Poincare 
recurrence theorem. 
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In this paper, we intend to give applied mathematics arguments for the validity of the 
famous Ergodic theorem in a class of polinomial non-linear and Lipschitz wave motions. 
Our approach is fully based on Hilbert Space methods previously used to study Dynam- 
ical System's in Classical Mechanics which by its turn simplifies enormously the task of 
constructing explicitly the associated invariant measure on certain Sobolev Spaces: the 
infinite - dimensional space of wave's motion initial conditions. This study is presented 
on the main section of this paper, namely section 4. 

In section 2, we give a very detailed proof of the RAGE'S theorem, a basic functional 
analysis rigorous method used in our proposed Hilbert Space generalization of the usual 
finite - dimensional Ergodic theorem and fully presented in section 3. 

In section 5, we complement our studies by analyzing the important case of wave - 
diffusion under random stirring. 

In appendices A and B, we additionally present important technical details related to 
the section 4 and they contain material as important as those presented in the previous 
sections of our work. 

2 On the detailed mathematical proof of the R. A.G.E. 
theorem. 

In this purely mathematical first section 2 of our study, we intend to present a detailed 
mathematical proof of the R.A.G.E theorem ([2]) used on the analytical proof of ours of 
the Ergodic theorem on section 3 for Hamiltonian systems of A-particles. 

Let us, thus, start our analysis by considering a self - adjoint operator £ on a Hilbert 
Space (H, (,)) where H C (C) denotes the associated continuity sub-space obtained from 
the spectral theorem applied to C We have the following result: 

Proposition (RAGE theorem) - Let ip e H C (C) and ip G H. We have the Ergodic 
limit m 



Proof - In order to show the validity of the above ergodic limit, let us re-write eq.(l) in 




(1) 
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terms of the spectral resolution of £, namely. 



T 



(■0, exp(itC)ip) 



dt 



L 



+00 

e""" d^E^) 

Here we have used the usual spectral representation of £ 

r+00 

(g,Ch)= / Xd x (g, Ej(X)h) 



x 



(2) 



(3) 

J —00 

with g E H and h G Dom (£). 

Let us remark that the function exp(i(A — fj)t) is majorized by the function 1 which, by 
its turn, is an integrable function on the domain [0, T] x R x R with the product measure 

dt 



T 



since 



T 



dt 



— ®dx(il>,E j (\)4>)®d ll (il>,E j ( t j,)4>) = (ip^y < 00 



T 



(4) 



(5) 



At this point, we can safely apply the Fubbini theorem for interchange the order of 
integration in relation to the t-variable which leads to the partial result below 

"+00 p+00 / i(\-n)T _ 



/+00 r+00 / 
-oo J —oo V 



d x {4>, Ej(X)ij) <g> d^, Ej(n)il>) 



i{\-n)T 

Let us consider two cases. Firstly we take tp — ip. In this case, we have the bound 



(6) 



e i{\-n)T _ l 



i{X-n)T 



< 



2 sen 



T) 



(A-/i)T 



< 1 



(7) 



If we restrict ourselvers to the case of A 7^ /i, a direct application of the Lebesgue 
convergence theorem on the limit T — > 00 yields that 1 = 0. 

On the other hand, in the case of A = \i we intend now to show that the set /i — X is 
a set of zero measure in R 2 in respect with the measure d\{ijj, Ej(X)ip) <S> d^ip, Ej(fj)ip). 
This can be seen by considering the real function on R 



/(A) = (^(A)V) 



(8) 
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We observe that this function is continuous and non-decreasing since ip G H c (£) and 
range Ej(Xi)C range Ej(X 2 ) for Ai < A 2 . Now, /(A) is an uniform continuous function on 
the whole real line R, since for a given e > 0, there is a constant A such that 



C0, Ej(-oo, -X)ip) < 



(i/>,Ej(-oo, A)V) < - 



(9) 



(10) 



Additionally, /(A) is uniform continuous on the closed interval [—A, A] and /(A) can not 
make variations greater than | on (— oo, — A] and [+A, oo), besides of being a monotonic 
function on R. These arguments show the uniform continuity of /(A) on whole line R. 
Hence we have that for a given e > 0, exists a 5 > such that 



(^,^(A')^)-(^^(A» 



< 



for 



A' - A" 

In particular for A' = A + 5 e A" = A — 5 



< 5 



(^,^(A + ^)-(^,^(A-^)< 



(11-A) 



(12) 



(13) 



As a consequence, we have the estimate 



-co /"A+<5 

d x {E 3 {X)^^) / d^EjfaW 

x. JX-S 



< 
< 



x ((^^.(A + W-^A-^)) 



x 



2 — 



< e 



(14) 



Let us note that for each e > 0, there is a set D$ = {(A,//) G -R 2 ; |A — //| < 5} 
wich contains the line A = ji and has measure less than e in relation to the measure 
d\(Ej(X)ijj,'ilj) <g> d fl (Ej(fj l )ip,tj;) as a result of eq.(14). This shows our claim that / = in 
our special case. 

In the general case of ip ^ ip, we remark that solely the orthogonal component on the 
continuity sub-space H C (C) has a non - vanishing inner product with exp(itC)ip. By using 
now the polarization formulae, we reduce this case to the first analyzed result of / = 0. 
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At this point we arrive at the complete R.A.G.E theorem ([2]) 



Theorem - Let K be a compact operator on (if, ( , )). We have thus the validity of the 
Ergodic limit 

i r T 

lim - / \\K exp{it£)ip\\ 2 H dt = (15) 

with ip G H c {£). 

We leave the details of the proof of this result for the reader, since any compact 
operator is the norm operator limit of finite - dimension operators and one only needs to 
show that 

lim )- I || V c„(e l S/>, e n )g n \\ H dt = (16) 
for c n constants and {e n }, {g n } a finite set of vector of (H, (, )) . 

3 On the Boltzman Ergodic Theorem in Classical 
Mechanics as a result of the R.A.G.E theorem 

One of the most important statement in Physics is the famous zeroth law of thermo- 
dynamics: "any system approaches an equilibrium state". In the classical mechanics 
frameworks, one begins with the formal elements of the theory. Namely, the phase-space 
R 6N associated to a system of iV-classical particles and the set of Hamilton equations 

3H . dH 

Pi = —5—; Qi = ^— ( 17 ) 

% dpi 

where H(q,p) is the energy function. 

The above cited thermodynamical equilibrium principle becomes the mathematical 
statement that for each compact support continuous functions C C (R 6N ), the famous er- 
godic limit should holds true ([3]). 

/ d 3N q(0)d 3N p(0) ( lim 1 f f(q(t);p(t))dt\ = V (f) (18) 

J RON [T^oo 1 J Q J 

where rj(f) is a linear functional on C c (R m ) given exactly by the Boltzman statistical 
weight and {q(t),p(t)} denotes the (global) solutions of the Hamilton equations (12). 



5 



We aim at this section point out a simple new mathematical argument of the fun- 
damental eq.(18) by means of Hilbert Space Techniques and the R.A.G.E'S theorem. 
Let us begin by introducing for each initial condition (g(0),p(0)), a function oj qo;po {t) = 
(q(t),p(t)) : here (q(t),p(t)) is the assumed global unique solution of eq.(17) with pre- 
scribed initials conditions. Let U t : L 2 (R m ) — > L 2 (R m ) be the unitary operator defined 
by 

(Utf)(q,p) = f(u qoPo (t)) (19) 
We have the following theorem (the Liouville's theorem) ([5]). 

Theorem 1. U t is a unitary one-parameters group whose infinitesimal generator is —iL, 
where —iL is the essential self-adjoint operator acting on C^°(R m ) defined by the Poisson 
bracket. 



3N 



(Lf)( P , q ) = {f,H}(a,p) = £ ( - %L?£) U, „) 



i=l 



dqi dpi dpi dqi 



(20) 



The basic result we are using to show the validity of the ergodic limit eq. (18) is the 
famous R.A.G.E's theorem exposed on section 1. 

Theorem 2. Let G H c (—iL)(L 2 (R 6N ) , here H c (—iL) is the continuity sub-space asso- 
ciated to self-adjoint operator —iL. For every vector f3 G L 2 (R 6N ), we have the result 

rT 



lim — 

T^oo T 

or equivalently for every ip G L 2 (R 6N ) 

rT 



dt = 0, 



lim — 



(J3,U t if>)dt={l3,F keii _ iL) if>) 



(21) 



(22) 



where P 



ker(-iL) 



is the projection operator on the (closed) sub-space ker(-iL). 



That eq.(22) is equivalent to eq.(21), is a simple consequence of the Schwartz inequality 
below written 



mu t ti>))dt 



< 



(((3,U^)) 2 d? 



ldt 



or 



- / ({3,U^)dt < - / ((P,U t ^)) 2 dt 



T 



(23) 



(24) 



As a consequence of eq.(23), we can see that the linear functional r)(f) of the Ergodic 
theorem is exactly given by (just consider {3(p, q) = 1 on supp of Pkcr(-iL)C0)) 

V(f)= [ dqdp¥ ker{ _ iL) (f)(q,p). (25) 

JR6N 

By the Riesz' s theorem applied to r](f), we can re-write (represent) eq.(25) by means 
of a (kernel)-function h v (H)(q,p), namely 

V(f)= f d m qd m p(f(p,q))h v{H) (p,q) (26) 

where the function h^H) (q, p) satisfies the relationship 



or equivalently h v (H){q,p) is a "smooth" function of the Hamiltonian function H(q,p), by 
imposing the additive Boltzman behavior for h v (H){q,p) namely, h v ^H 1 +H 2 ) — ^(Hi) - h v (H 2 )j 
one obtains the famous Boltzman weight as the (unique) mathematical output associated 
to the Ergodic Theorem on Classical Statistical Mechanics in the presence of a thermal 
reservoir ([4]). 

h v{H) (q,p)=eM-PH(q,p)}/ J d 3N qd m p exp{-/3if(g,p)} (28) 

with (3 a (positive) constant which is identified with the inverse macroscopic temperature 
of the combined system after evaluating the system internal energy in the equilibrium 
state. Note that ||/wh) ||l 2 — 1 since ||^(/)|| = 1. 

A last remark should be made related to eq.(28). In order to obtain this result one 
should consider the non-zero value in ergodic limit 

1 f T 

Inn - / dt(U t ij)(q,p)h v{H) (q,p) = {h v(H) ,F 

ker (—iL) W) (29) 

or by a pointwise argument (for every t) 

ker (— iL)i 

(30) 

that is 

h v(H) e¥ kcr{ _ lL) ^ Lh = 0. (31) 
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4 On the invariant Ergodic functional measure for 
some non- linear wave equations 



Let us start this section by considering the discreticized (TV-particle) wave motion Hamil- 
tonian associated to a non-linear polinomial wave equation related to the vibration of a 
one-dimensional string of lenght 2a 

N 



2k 



(32) 



with the initial one boundary Dirichlet conditions 



?i(0) = a 

Pi(0) =Pi 

qi(-a) = qi{a) = (33) 

Here k denotes a positive integer associated to the non-linearity power and g > the 
non-linearity positive coupling constant. 

Since one can argue that the above Hamiltonian dynamical system posseses unique 
global solutions on the time interval t G [0, 00) one can, thus, straightforwardly write 
the associated invariant (Ergodic) measure on basis of the Ergodic theorem exposed in 
section 2 restricted to the configuration space after integrating out the term involving the 
canonical momenta 

lim 1 f dt F( qi (t), • • • , q N (t)) = [ F( Ql , • • • , q N )S inv ^(qi, •••,<?*) (34) 

T^oo 1 J J RN 

here the explicitly expression for the Hamiltonian system invariant measure in configura- 
tion space is given by [with (3 = 1] 

4 mV) ■ • • ,Qn) = exp I — ^ (77) (fe+i -9i)Va 2 ) + J^i) 2 * J x (dqi ■ ■ ■ dq N ) 

(35-a) 



Z<°> = / d^v( qi ,--- ,q N ) (35-b) 
Jr n 

Now it is a non-trivial result on the theory of integration on infinite dimensional 
spaces that the cilindrical measures (normalized to the unity !) converges in a weak 
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(star) topology sense to the well - defined non -linear polinomial Wiener measure at the 

N L 

continuum limit N — > oo, or ^() — > J_ L dx() and qi(t) — > U(x,t) (see appendix B for the 

1=1 

relevant mathematical arguments supporting the mathematical existence of such limit ). 

As a consequence of eqs.(35) at the continum limit, we have the infinite - dimen- 
sional analogous of Ergodic result for the non-linear polinomial wave equation (with the 
Boltzman constant equal to unity - see eq.(28) - section 3 of this study) 

1 f T 

lim- / F(U(x,t))dt 

= [ d Wtener fi[U {0 \x)]x F(U {0 \x))exp{-^ I" dx(U {0) (x)) 2k \ (36) 

JC a (-a,a) I 2 ^ J -a J 

where the wave field U(x ) t) E L°°([0, oo), L 2k ~ l ((— a, a))) satisfies in a generalized weak 
sense the non-linear wave equation below (see appendix A). 

^ = «^> +s(t , (M)) ~ (37) 

U(x,0) = U<®(x)eH 1 ((-a,a)) (38) 

U t (x,0) = (39) 

U(-a,t) = U(a,t) = (40) 

d Wiener n[U^(x)} denotes the Wiener measure over the uni-dimensional Brownian field 
end-points trajectories {U^{x)\ with U<®(-a) = U^(a) = 0} and F(U^(x)) denotes 
any Wiener - integrable functional. 

Note that C a (—a,a) denotes the Banach space of the a-Holder continuous function 
(for a <\) which is obviously the support of the Wiener measure d Wiener fi[U^ (x)]. 

It is worth remark the normalization factor Z of the non-linear Wiener measure im- 
plicitly used on eq.(36). 

± / d w ~^\x)} exp {-J- f dx{uV\x))A = 1 (41) 

Let us pass to the important case of existence of a dissipation term on the wave 
equation problem eq.(37) - eq.(40) with v the positive viscosity parameter and leading 
straightforwardly to solutions on the whole time propagation interval [0, oo) 



U(x,0) = U<®(x) eH\(-a,a)) (43) 

U t (x,0) = U {1) (x) e L 2 ((-a,a)) (44) 

It is a simple observation that there is a bijective correspondence between the solution 
of eq.(42) with the Klein - Gordon like wave equation for the re-scaled wave field (3(x, t) = 
e^U (x, t), namely 

d 2 /3(x,t) d 2 (3(x,t) fv 2 \ n . . (2fc-i)w , 

1 ^(x,t)+^e — —((3(x,t)) 2k 1 (45) 



9 2 t a% v 4 

/?( x ,0) = f/ (0) (a:) =^((-0,0)) (46) 

(3 t (x,0) = U^(x)-^°\x)=0 (47) 

/3(-a,t) = /3(a,*) = (48) 

As a result we have the analogous of the Ergodic theorem applied to eq.(45) and a 
sort of Caldirola-Kanai action is obtained ([5]), a new result of ours in this non-linear 
dissipative case (with the path-integral identification x — * a; U^(x) = X(a)) 



hmi f dt F(e + ^U(x,t)) 

= [ i d Wiener n[X{a)\ F(X(a)) 
xexp{-| £^ e -^x(X(a))-} 



xexp \ + W j a d ^ x ^)) 2 ) ( 49 ) 

Concerning the higher-dimensional case, let us consider A a strongly positive elliptic 

operator of order 2m associated to the free vibration of a domain Q in a general space 
R u , and satisfying the Garding coerciviness condition 

A= £(-l) |a| ^(a Q/3 (x))L>f (50) 

|a|<m 
\P\<m 

with the Sobolev spaces operator's domain 

D(A) = H 2m (Q) n H™ (Q). (51) 
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Note that exists a constant C (Q) > (C (Q) — > if vol (Q) — > oo) such that the Garding 
coerciviness condition holds true on D(A) 

Real (AU,U) LHn) > C (Q)\\f\\ H 2 m{n) (52) 

We can thus associate a Lipschitz non-linear external vibration on the domain Q as gov- 
erned by the wave equation below. 

^^ = (AU)(x,t) + G'(U(x,t)) (53) 

U(x, 0) = f/ (0) (x) e H 2rn (O) n H™{Q) (54) 
U t (x,0) = g(x) e L 2 (n) (55) 

Here G(z) denotes an differentiable Lipschitizian function, like G(z) = A sen (Bz) , e~ 7Z , 
etc... and thus, leading to the uniqueness and existence of global weak solutions on 
C([0, oo), H 2m (tt) n H™(tt)) by the Fixed point theorem ([6]). 

Analogously to the discretization j A^-particle technique of the 1 + 1 case, one obtains 
as a mathematical candidate for the Ergodic invariante measure, the following rigorously 
measure on the Hilbert Space of functions H 2m (Vt) (with (3=1) 

d {inv) i2((f)(x)) = d A/ 2((j)(x)) x e -3 fn dVx G ^(*)) (56) 

Here dAfi(4>(x)) denotes the Gaussian measure generated by the quadratic form asso- 
ciated to the operator A on its domain D(A) for real initial conditions functions <f>(x). 

[ d A /x(0(x))(0(x 1 )0(x 2 )) = (A~ 1 )(x 1 ,x 2 ) (57) 

Note that the Green function of the operator A belongs to the trace class operators 
on H 2m (Q), which results on eq.(57) through the application of the Minlo's theorem ([7]). 

At this point, we remark that di nv /j,((f)(x)) defines a mathematical rigorous Radon 
measure in H 2m (Q) since the function (functional) 

F(<f>(x))=exp\~ f dTx G(<f>(x))\ (58) 



2__ 

belongs to L 1 (f/" 2m (J7), dAfJ>(4>(x))) due to the Lipschitz property of the non-linearity 
G(z), i.e. there is constants c + and C~, such that 

c" cf){x) < G(<j>(x)) < c + (f)(x) (59) 
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As a consequence, we have the obvious estimate below 

e -c+f n d"uxcf>(x) <e ~U n <Tx G(4>(x)) < g-c" /„ d"x^x) 



(60) 



so by the Lebesgue comparasion theorem, we get the functional integrability of the non- 
linearity term of the ergodic invariant measure eq.(56): 



>H 2m (VL) 



d A fi(<j)(x))e z J » 



i f n d»x G(4>(x)) 



< e +|(^+c + ) 2 Tr fl2m(n) (A- 1 ) < 



OO 



(61-a) 



and leading thus to the mathematical validity of eq.(56). 

At this point it is worth remark that one can easily generalize the type of our allowed 
Lipschitz non-linearity for those of the following form 

/ d v x G(<j>(x))\ > (+ / d u x^-(<P) 2 (x)+ [ d v xC{x)<f>{x)\ (61-b) 
Jn J I Jn 2 J n J 

with 7(x) and C(x) real positive functions. 

This claim is a result of the straightforward Gaussian - infinite - dimensional integra- 
tion below 



H 2m (n) 



^(0)exp<j- / drx^-(H>){x) + I <FC{x)<t> {x) 



det 2 
1 



< 



exp < + 



< oo 



n * Jn 
1 + J d v yA-\x-y)- 



i{y) 



d v x / d»yC{x)-A-\x,y)C{y) 



n 



(61-c) 



In the usual (Euclidean) Feynman path-integral formal notation, the Ergodic theorem 
takes the physicist's form after re-introducing the reservoir temperature parameter j3 = 
1/kT. 

lim i [ F(U(x,t)) = \\ D F [ V (x)}e'^^ dVx ^ {x) - {Av{x))) x F(<p(x)) 
with the normalization factor 



Z = 



ff 2m (Q) 



D F [^(x)\e^ 13 ^ {x){Aip){ - x))dVx e~ p ^ d " x GMx)) 



(63) 
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The functionals F((p(x)) are objects belonging to the space L 1 (if 2m (fi), d A ((p(x))) C 
C{H 2m {Q),R). 

At this point let us remark that by a direct application of the mean value theorem 
there is a sequence of growing times {t n }„ e z (with t n — > oo) such that 

lim F(0(x,t n )) = - / L» F [^(x)]e-^io^(^)W 

x F{ V {x))e-P^ dVxG ^ x)) (64) 

It is worth recall that the (positive) parameter (5 is not determined directly from the 
parameters of the mechanical underlyng system and can be thought as a remnant of 
the initial conditions "averaged out" by the Ergodic integral (see section 3) and some- 
times related to the presence of a intrinsic dissipation mechanism involved on the Ergodic 
- time average process responsible for the extensivity of the system's thermodynamics 
(see eq.(28) - section 3). For instance, at the limit of vanishing temperature (absolute 
Newtonian-Boltzman-Maxwell zero temperatire (3 — > oo and somewhat related to the 
so called ergodic mixing microconomical - Ensemble ), on can see the appearance of an 
"atractor" on the space H 2m (Q) fl H™(Q) and exactly given by the stationary strong 
unique solution of the wave equation eq.(53) since, we always have an soluble Elliptic 
problem on H 2m (Q) fl H™(Q) by means of the Lax-Milgram-Gelfand theorem ([6]). 

(¥)W=G'(/W) (65) 

if*(x) = (66) 
an 

note that the stationary strong solution ip*(x) is not a real atractor in the sense of the 
treory of Dynamical systems, since we can only say that there is a specific sequence of 
growing {t n } ne z, with t n — > oo, such that at the leading asymptotic limit (5 — > oo we have 
the time-infinite limit for every solution of eqs.(65) - (66). 

lim F(<p(x,t n )) = F(<p*(x)) (67) 

A (somewhat) formal path-integral calculation of the next term on the result eq.(41) 
for F(z) being smoth real- valued functions can be implemented by means of the usual 
path-integrals Saddle-point techniques applied to the (normalized) functional integral 



13 



eq.(62) ([8]). 



lim F(<p(x,t n )) = F(<p*(x)) 



o=0 



F(<p*(x)) + \ \ ^V(x)) x Tr H2m(n) 



(68) 



/3 

{5 2 F 



, i ^ 2 C/ * 



+ 0(/T 2 ) 



(69) 

Another important remark to be made is related to the existence of time-independent 
simple constraints on the discreticized Hamiltonian eq.(32), namely 



H {a) (q ir -- ,q N ) = 



a — 1, • • • , n 



(70) 



In this case, one should consider that the continuum version of eq.(70) leads to 
(strongly) continuos functionals on H 2m (Vl) fl H™(Vt) 



# (a) (0OM)) = 



a — 1, • • • , n 



(71) 



They have a direct consequence of restricting the wave motion to a closed non-linear 
manifold of the original vector Hilbert Space H 2m (Q) fl H™(Tl). 

Now it is straightforward to apply the Ergodic theorem on the effective mechanical 
configuration space and obtain, thus, as the Ergodic invariant measure the constraint 
path-integral. 

i r T 

lim-/ dt F(<f>(x,t)) 

= — f D F [<p(x)]e-&fn' Px M x) ' Av ( x » e -PSnd v xG(<P{x)) 

Z JH' 2m (n) 



x 



(72) 



va=l 



where 5^ F \-) denotes the usual formal delta - functionals. 
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5 An Ergodic theorem in Banach Spaces and Appli- 
cations to Stochastic-Langevin Dynamical Systems 

In this complementary section 5, we intend to present our approach to study long-time 
/ ergodic behavior of infinite-dimensional dynamical systems by analyzing the somewhat 
formal diffusion equation with polinomial terms and driven by a white noise stirring ([4]). 

In order to implement such studies, let us present the author's generalization of the 
R.A.G.E theorem for a contraction self-adjoint semi-group T(t) on a Banach space X. 
We have, thus, the following theorem. 



Theorem 3. Let f £ X . We have the ergodic generalized theorem 

m — / 

T- 



lim i f dt(e- tA f) = F ker(A) (f) (73) 



where A is the infinitesimal generator ofT(t). 

Let us scketch the proof of the above claimed theorem of ours. 

As a first step, one should consider eq.(73) re- written in terms of the "resolvent opera- 
tor of A" by means of a Laplace Transform (The Hile-Yosida - Dunford Spectral Calculus) 

^J[ dt {^L dze " ((z+Arlf) }- (T4) 

Now it is straightforward to apply the Fubbini theorem to exchange the order of 
integrations (dt,dz) in eq.(74) and get, thus, the result 

lim 1 f dt(e- tA f) = lim +((* + A)^f). (75) 

The z — > CT limit of the integral of eq.(75) (since Real (z) C p(A) C (— oo, 0) ) can be 
evaluated by means of saddle point techniques applied to Laplace's transforms. We have 
the following result 



limKz + A)- 1 /) 

poo 

= lim / e- zt (e~ tA f) dt 

*-+ 0- Jo 

= lim(e- M /)=P M ,)(/). (76) 
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Let us apply the above theorem to the Lengevin Equation. Let us consider the Fokker- 
Planck equation associated to the following Langevin Equation 

f(') = -f WW) + ,'(*) (77) 

where {rf(t)} denotes a white- noise stochastic time process representing the thermal cou- 
pling of our mechanical system with a thermal reservoir at temperature T. Its two point 
function is given by the "Fluctuation-Dissipation" theorem 

( V i (t)TfV)) = kT6(t-lf). (78) 

The associated Fokker-Planck equation associated do eq.(74) has the following explic- 
itly form 

^(q\ <?; t) = +(kT)A qi P(q\ q\ t) + V w (V ft V • P(q\ (f; t)) (79) 

kP^^^f'^-?-). (80) 

By noting that we can associate a contractive semi-group to the initial value problem 
eq.(80) in the Banach Space L 1 (_R 3Ar ), namely: 

J P{q\q\t)f(<f)dq l = (e- tA f). (81) 

Here the closed positive accretive operator A is given explicitly by 

-A(.) = *rA fli (.) + V ft [(V ft V)(.)] (82) 

and acts firstly on X = Coo(R 3N ). It is instructive to point out that the perturbation 
accretive operator B = V qi ■ [(V gi V) ■ (•)] on ^(i? 3 ^) with V(q l ) E C c °°(i? 3JV ) is such 
that it satisfies the estimate on S(R 3N ) : \\B fWs^N) < a\\ Af\\ S (R3N) + 6||/||s(.r3.iv) with 
a > 1 and for some b. As a consequence A(-) = kTA qi (-) + (B(-)) generates a contractive 
semi-group on Coo(-R 37V ) or by an extension argument on the whole L 1 (i? 37V ) since the 
L 1 -closure of Coo(R 3N ) is the Banach space L 1 (i? 3Ar ). 

At this point we may apply our Theorem 3 to obtain the Langevin-Brownian Ergodic 
theorem applied to our Fokker-Planck equation 



hm 1 f dt( [ dq(e- tA f)(q)) = [ dqF ker{A) (f)(q) = [ d N <f f {(?)^ {(?) 

1 J0 \JR 3N J JB? N JR N 

(83) 
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where the equilibrium probability distribution is given explicitly by the unique normal- 
izable element of the closed sub-space ker(A). 

O = +(kT)A q F eq (q l ) + Vf[(VfV)Ftf)] (84) 

or exactly, we have the Boltzman's weight for our equilibrium Langevin-Brownian prob- 
ability distribution 

(85) 



P e W=exp j- — (V(<?)) 

For the general Langevin equation in the complete phase space {qi,Pi} as in the bulk 
of this note, one should re-obtains the complete Boltzman statistical weight as the equi- 
librium ergodic probability distribution. 



with the normalization factor 



df 



R N 



R N 



(86) 



(87) 



Let us apply the above exposed result for the following non-linear (polinomial) stochas- 
tic diffusion equation on a one-dimensional domain Q = (—a, a) with all positive coeffi- 
cients {Xj} an the non-linear polinomial term 

dU(x,t) ld 2 U(x,t) 2/, '" J 



dt 



+ - 



d 2 x 



"(E ^U j )(x,t)+ V (x,t) 



(88) 



3=0 
j=odd 







(89) 
(90) 



U(x,t) 

an 

U(x,0) = f(x) G L 2 (9). 

Here {r)(x,t)} are samplings of a white-noise stirring. 

In this case, one can show by standard techniques that the global solution U (x, t) G 
C([0, oo), L 2 {— a, a)) and by using the discreticized approach of section 4, the invariant 
Ergodic measure associated to the non-linear diffusion equation is given by the mathe- 
matically functional equilibrium - Langevin Brownian probability distribution 

/ 2 k . \ 



P e i(U{x)) = exp 



dx 



E j ui < 

3=0 J 

j even 



X) 



J 



x d_ 1 d2 n(U{x)) (91) 

2 
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Here d_ L d 2 (i(U(x)) denotes the Gaussian functional measure (normalized to unity) on 



2 dx^ 



the Sobolev Space if 1 (— a, a) space. In the Feynman path-integral notation (with (3=1) 

1 f T 

lim - / dtF(U(x,t)) 



L 



D F [U(x)}e-y-«^ 2 W 



j=0 J 

x e j=even x F[U(x)] (92) 

Note that eq.(91) defines a rigorous mathematical object since it does not need to be 
renormalized from a calculational point of view on the Sobolev space Hi(—a,a) ([10]) as 
its non-linear exponential term is bounded by the unity (Lebesgue theorem). It is worth 
point out that for non-linearities of Lipschitz type on general domains QCR U , one can see 
that the support of the Gaussian measure is readly H 2 (Q) and the equilibrium Langevin 
- Brownian probability distribution 



P e «(U(x)) = exp j-^ d D x G{U(x))^ x d { _i A) v(U(x)) 



(93) 



is a perfect well-defined Randon measure on H 2 (Vl) by the same Lebesgue convergence 
theorem. 

6 The existence and uniqueness results for some non- 
linear wave motions in 2D 

In this technical section 6, we give an argument for the global existence and uniqueness 
solution of the Hamiltonian motion equations associated firstly to eq.(32) - section 3 and 
by secondly to eq.(53) - section 3. Related to the two-dimensional case, let us equivalently 
show the weak existence and uniqueness of the associated continuum non-linear polinomial 
wave equation in the domain (—a, a) x R + . 

d 2 U(x,t) d 2 U(x,t) , IJTI ,. 2fc _ 1 



d 2 t d 2 x 



+ g(U(x,t)Y k - 1 = (94) 
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U(-a,t) = U(a,t) = 
U(x, 0) = U (x) G H 1 ( [-a, a] ) (95) 

U t (x,0) = U 1 (x)EL 2 ([-a,a]) (96) 
Let us consider the Galerkin approximants functions to eq.(95) - eq.(96) as given below 

U n (t) = J2 Ui(t) sen f-x) (97) 
t=i \ a J 

Since there is a 70(a) positive such that 

d 2 _ " 



U n ,U n ) > lo(a) (U n ,U n )m([-a,a}) (9£ 

L2([-a,a]) 



we have the a priori estimate for any t 

< y?(t) < (p(0) (99) 

with 

^(0 = ^lFnWlli 2 +7o(a)||C7 n ||^ 

+ ^\\Un\\% (100) 

As a consequence of the bound eq.(100), we get the bounds for any given T (with Ai 
constants) 

sup ess||I7 n ||^i ( _ } < Ax (101) 

0<t<T 

sup ess\\U n \\ 2 L 2, a) < A 2 (102) 

0<t<T 

sup ess\\U n \\% ( _ aa] < A 3 (103) 

0<t<T 

By usual functional - analytical theorems on weak-compactness on Banach-Hilbert 
Spaces, one obtains that there is a sub-sequence U n (t) such that for any finite T 



U n (t)(weak-star)U(t) in L°°([0, T], H\-a, a)) (104) 
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U n {t) (weak - star) v(t) in L°°([0, T], L 2 {-a, a)) (105) 



U n (t)(weak- star)p(t) in L°°([0, T], L 2fc (— a, a)) (106) 

At this point we observe that for any p > 1 (with constants) and T < oo we have 
the relationship below 

f \\Un\\ P m{ _ a , a) dt < T{Arf # f \\U n \\l H _ aA) < M (107) 

J J 

/ T Pnll^^)* < T(A 2 )i ^ ^ lFHI 2fc( _ a , a) ^ < M (108) 
since we have the continuous injection below 

H\-a,a) L 2 (-a,a) (109) 

L 2fe (-a,a) L 2 (-a,a) (110) 

As a consequence of the Aubin-Lion theorem ([6]), one obtains straightforwardly the 
strong convergence on L p ((0,T), L 2 (—a, a)) togheter with the almost everywhere point 
wise equalite among the solutions candidate 

U n — >U(t)=p(t)= [ v{s)ds (111) 

Jo 

By the Holder inequalite applied to the pair (q, k) 

\\U n - U\\ L , { - a , a) < Wn-UWvUa) 1 1 U * ~ ^ I ^ 

with < e < 1 

1 1-9 9 , x 

r~ + Tk < n3 > 

in particular with q — 2k — 1, one obtains the strong convergence of U n (t) in the general 
Banach space L°°((0, T), L 2k -\-a, a)), with = ^ (|gf ) 

As a consequence of the above obtained results, one can pass safely the weak limit on 

C™((0,T),L 2 (-a,a)) 

( dp \ 

= ^(U,v)^-a,a) + \-^P' V ) L2 _ +9(U 2k ~\v) LH _ a , a) = (114) 
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for any v G C°°((0, T), L 2 (-a, a)) 

At this point we sketchy a somewhat rigorous argument to prove the problem's unique- 
ness. 

Let us consider the hypothesis that the finite function 

w= gfeerM („ 5) 

(E/"(a;,t) -v(x,t)) 

is essentially bounded on the domain [0, oo) x (—a, a) where U(x,t) and t) denotes, 
two hypothesized different solutions for the 2D-polinomial wave equation eq.(94) - eq.(96). 
It is straightforward to see that its difference W(x,t) = (U — v)(x, t) satisfies the "Linear" 
wave equation problem 

d 2 W d 2 W 

___ + (oH0(M) = (116) 
W(0) = (117) 
W t (0) = (118) 
At this point we observe the estimate (where H 1 <^-> L 2 \) 

II II (\\T*r\\ \\ dW \\ 

< \\a\\L<x>((0,T)x(-a,a)) x I || W || L 2 (-a,a) x \\ —jjT \\ L 2 )-a,a) 

^ M (\\^\\lH-a,a) + \\W\\h { -a,a)^ 

<M(||^||i 2( _ a , a) + ||^||^i ( _ a , a) ) (119) 
which after a application of the Gronwall's inequalite give us that 

J t W WU-a,a) + \\W\\ 2 m) ^ (t) < (||^||! 2( _ a , a) (0) + ||W||^(_ a , a) (0)) = (120) 

which proves the problem's uniqueness under the not proved yet hypothesis that in the 
two-dimensional case (at least for compact support infinite different iable initial conditions) 

sup \\a(x,t)\\<M (121) 

x£(— a, a) 

*e[o,oo) 
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As a last comment on the 1 + 1-polinomial non-linear wave motion, we call the reader 
attention that one can easily extend the technical results analyzed here to the complete 
polinomial non-linearity 

x2jt a2TT i <2k 

Cj U 2k - j = C 1 U 2k - 1 + C 3 U 2k ~ 3 + --- (122) 



d 2 t d 2 x 

3=1 
j,odd 



with the set of couplings {Cj} belonging to a positive set of real number. In this case we 

have that the general solution belongs a priori to the space of functions f] L°°((0,T), 

3 i 

j<2k 

L 2k -i{-a,a)). 



22 



Appendix A 

On sequences of Randon measures on Functional Spaces. 



Let us consider a completely regular topological space X and the Banach space of 
continuous function space there 

C(X) = {/ : X — > R, f continuos on the X- topology with ||/||oo = sup \f(x)\ < 00} 

xex 

(A-l) 

Associated to the Banach Space C(X), there is its dual M(X) which is the set of 
Random measures on X with norm given by its total variation << ji »= \fJ>\(X) (where 
l*eM(X)). 

It is easy to see that if f3(X) is the Stone-chech compactification of X, we have the 
isometric immersion Ai(X) <^-> /A(/3(X)) 

We have, thus, the important theorem of Prokhorov on the accumulation points of 
sets of Random measures ([7]). 

Theorem (Prokhorov) - If a family of measures on M(X) is such that it satisfies 
the following conditions 

a) Bounded uniformely 

«H a »=\n a \{X)<M (A-2) 

c) Uniformely regular - For any given e > 0, there is a compact set C X such that 

f(x) dfx a (x) 
for any f(x) G C(X) satisfying the condition 

(supp/(*))ntf (e) = {0} (A-3) 

We have that the set of Random measures is a relatively compact set on the weak-star 
topoly on A4(X). Namely, for any f(x)C(X), we have that there is a measure v e A4(X) 
such that for a given sequence {at} 



x 



^ £|| J I ! A. 
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lim sup / f(x)dfj, aic (x) = / f(x)du(x) (A-4) 

{a K } JX JX 

In order to apply such theorem to the set of Random measures as given by eq.(35), 
we introduce the completely Regular space formed by the one-point compactification of 
the Real line 



X = U.(R)i (orX = l[j3(R)) (A-5) 

and note the strict positivity of the integrand (the path- integral wheight) on eq.(35) 
and its boundedness by the Gaussian free term since exp | — |^f^ Zj(<?«) 2fe | < 1 f° r 

N 

{?.}en(4 

1=1 

We note either the chain property of the invariant measure set ensuring us the measure 
uniqueness limit on eq.(B-4) when applied to eq.(35) - section 4, namely: 

/ f( Ql , • • • , gjv -i) d^Vfoi, ■■■,Qn)> [ f(qi, • • • , q N -i) d^^(q u • • • , q N ^) 
Jr n Jr 1 *- 1 

(A-6) 

Finally, we remark that the well-defined non-linear Wiener functional measure coin- 
cides with those finite-dimensional as given by eq.(35-a) - eq.(35-b) when it is restricts 
to the set of continuous poligonal paths and, thus, ensuring the convergence of the mea- 
sure set eq.(35-a) - eq.(35-b) to the full Wiener path measure eq.(36), since its functional 
support is contained on the set of continuous functions on the interval [—a, a]. 
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